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■ We study the dynamics of one-particle and few-particle billiard systems in containers of var- 
ious shapes. In few-particle systems, the particles collide elastically both against the boundary 
and against each other. In the one-particle case, we investigate the formation and destruction 
of resonance islands in (generalized) mushroom billiards, which are a recently discovered class of 
Hamiltonian systems with mixed regular-chaotic dynamics. In the few-particle case, we compare the 
dynamics in container geometries whose counterpart one-particle billiards are integrable, chaotic, 
and mixed. One of our findings is that two-, three-, and four-particle billiards confined to contain- 
ers with integrable one-particle counterparts inherit some integrals of motion and exhibit a regular 

\ partition of phase space into ergodic components of positive measure. Therefore, the shape of a 

, container matters not only for noninteracting particles but also for interacting particles. 

^ ; PACS numbers: 05.45.-a, 05.10.-a 

In this paper, we conduct a numerical investigation of one-particle systems (billiards) with regular, 
(S| ' chaotic, and mixed (regular-chaotic) dynamics and of small numbers (two, three, and four) of elasti- 
^ . cally colliding particles (balls) confined to the same billiard tables. Thus, this report differs essentially 

■ from traditional numerical studies of systems of interacting particles, which deal with many particles 
confined to containers of the simplest shapes (usually boxes or boxes with periodic boundary condi- 

QQ ■ tions [i.e., tori]). Using the simple example of hard balls in a circle, we demonstrate that systems of 

. interacting particles need not be ergodic. Therefore, the recently-discovered typical inhomogeneity 

■ of stationary distributions of noninteracting particles in containers may well be in effect for interact- 
. ing particles as well. We also compare the dynamics of few-particle systems with their one-particle 

counterparts when the dynamics of the latter are regular, chaotic, and mixed. 



INTRODUCTION 



> 

\ Two major 20th century discoveries completely transformed scientists' understanding of nonlinear phenomena. One 
^ ■ was Kolmogorov-Arnold-Moser (KAM) theory, which demonstrated the stability of regular dynamics for small pertur- 
bations of Hamiltonian systems |5|, l20|, • The other was the theory of stochasticity of dynamical systems (loosely 
called "chaos theory"), which demonstrated the stability of strongly irregular dynamics under small perturbations 

Typical Hamiltonian systems have mixed dynamics, with islands of stability ("KAM islands") situated in a chaotic 
sea. However, rigorous mathematical investigations of such systems are notoriously difficult because different analytical 
methods have been developed for systems with fully regular or fully chaotic dynamics. Both approaches fail at the 
boundaries between chaotic and regular regions. Numerical investigations of systems with mixed dynamics are also 
difficult for the same reason (small islands are not easy to find/observe numerically). 

There have been numerous attempts to find Hamiltonian systems with mixed dynamics (divided phase space) that 
allow an exact, rigorous analysis 0, O, Recently, unexpectedly simple and visual examples of such systems 

were found in the form of mushroom billiards, whose geometry generalizes the long-studied stadium billiard 
From a mathematical perspective, the discovery of mushroom billiards has now made it possible to address some 
delicate questions about the dynamics of systems with coexisting KAM islands and chaotic regions 0, "o] . Such 
theoretical studies can be readily tested experimentally in many physical situations. For example, two-dimensional 
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billiards in essentially arbitrary geometries corresponding; to systems with integrable, chaotic, and mixed dynamics 
can be constructed using microwave cavities I^ISH^, quantum dots |2^], and atom optics ^,0,|3^. The problems 
that can be studied by experimentalists using billiard systems are both fundamental and diverse, ranging from the 
decay of quantum correlations 3] to investigations of the dynamics of Bogoliubov waves for Bose-Einstein condensates 
confined in various billiard geometries 36]. 

In this work, we investigate chaos-chaos transitions in classical billiard systems. To do this for a given billiard, 
we gradually perturb its geometry, just as one would for the order-chaos transitions that have long been studied in 
Hamiltonian systems. In this case, however, the billiard's phase space has chaotic regions both before and after the 
perturbation rather than just after it. We examine the proliferation and destruction of KAM islands that result from 
these perturbations and consider the concomitant disappearance and appearance of chaotic regions in phase space. 
In particular, we investigate the example of one-particle generalized mushroom billiards. We also study few-particle 
billiards, in which the confined particles collide elastically not only against the boundary but also against each other. 
Towards this end, we study them in geometries whose associated one-particle billiard dynamics are integrable, chaotic, 
and mixed. We consider particles of different sizes to examine the limit in which the particles become smaller and 
the few-particle system becomes more like the associated one-particle billiard. 

Traditional mathematical studies of systems of particles deal either with just one particle (billiards and pertur- 
bations/modifications thereof) or with an infinite number of particles. Numerical studies have been concerned with 
hundreds or thousands (or more) interacting particles Q, 12^ . Another tradition for studies of systems of inter- 
acting particles is to focus on the potential governing interactions between particles rather than on the shape of the 
containers in which the particles are confined; systems of interacting particles have traditionally been considered in a 
box or a torus. 

However, from both the theoretical and practical standpoints, systems of just a few particles are getting increasing 
attention, in large part because of nanoscience and possible applications in nanotechnolo gy. From an experimental 
perspective, collisions between particles can be studied using the framework of cold atoms |l5l On the theoretical 
side, it was recently demonstrated that stationary distributions of systems of noninteracting particles in containers 
with nontrivial shapes are typically very inhomogeneous 9]. A natural question that consequently arises is whether 
such phenomena may occur in systems of interacting particles. In this paper, we give a positive answer to this question 
by demonstrating that the same effect is present in systems of two, three, and four hard balls in a circle. 

The structure of this paper is as follows. First, we examine the dynamics of one-particle billiards. We give 
some background information on integrable and chaotic billiards before turning our attention to billiard systems 
with mixed regular-chaotic dynamics. We consider, in particular, billiard geometries shaped like mushrooms and 
generalized mushrooms. We then study the dynamics of two-particle billiards in containers whose corresponding one- 
particle billiard dynamics are integrable, chaotic, and mixed. We subsequently examine three-particle and four-particle 
billiards. 



ONE-PARTICLE BILLIARDS 

Classical (one-particle) billiard systems are among the best-studied Hamiltonian dynamical systems US EH. The 
dynamics of a billiard are generated by the free motion of a particle (usually taken to be a point) inside a closed 
domain Q with piecewise smooth boundary dQ in Euclidean space R^. (Most studies concentrate on the case d = 2.) 
The confined particle collides elastically against the boundary, so the angle of incidence equals the angle of refiection. 

The billiard fiow has phase space Ad = {{q^v) : q G \\v\\ = 1}. The unit normal vector (pointing towards the 
interior of the billiard) to dQ at the regular point q G dQ is denoted n{q). Also, let z := {q^ v) denote a point in phase 
space. 

A natural projection of M. into its boundary is M = {{q^v) : q G dQ^ \\v\\ = 1, {v^n{q)) > 0}, where (•, •) denotes 
the standard inner product. Define the billiard map T by T{q^v) = (g'i,vi), where qi is the point of dQ at which 
the oriented line through {q^v) first hits dQ and vi = v — 2{n{qi)^v)n{qi) is the velocity vector after the refiection 
against the boundary at the point qi. (The billiard map is not defined if qi is at a singular point of dQ.) When d = 2, 
this map is parametrized according to the arclength s along dQ (measured from an arbitrary point on dQ) and the 
angle (p between v and n{q). Hence, cp G [— 7r/2,7r/2] and {n{q)^v) = cos(p. The billiard fiow conserves the volume 
measure dv m. M. The corresponding invariant measure for T is d/j. = A cos (pdsdip, where A = l/[2 length(9(5)] is a 
normalization constant. 
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Billiards with Fully Integrable or Chaotic Dynamics 

A curve F C Q is called a caustic of the billiard if whenever any link of some trajectory is tangent to F, then all 
other links of the same trajectory are also tangent to F. It is well-known that the configuration space of a billiard 
in a circle is continuously foliated by circles concentric to dQ. A circular billiard is thus completely integrable. 
Elliptical billiards are similarly integrable, although they have two continuous families of caustics instead of just one. 
One family is formed by trajectories tangent to confocal ellipses, and the other is formed by trajectories tangent to 
confocal hyperbolae. In the limit of zero eccentricity, the "hyperbolic" caustics disappear and the elliptical caustics 
become concentric circular caustics. Semicircles and semiellipses are likewise integrable 19j. 

Examples of chaotic billiards include dispersing billiards such as the Sinai and diamond containers and billiards 
with focusing boundaries such as stadia. Neighboring parallel orbits diverge when they collide with dispersing compo- 
nents of a billiard's boundary. In chaotic focusing billiards, on the other hand, neighboring parallel orbits converge at 
first, but divergence prevails over convergence on average. (Hence, billiards with focusing boundaries may be viewed 
as occupying an intermediate position between dispersing and integrable billiards.) Divergence and convergence are 
balanced in integrable billiards. 

The investigation of chaotic billiards has a long history. The fiow in classical chaotic billiards is hyperbolic, ergodic, 
mixing, and Bernoulli 1^3 • Autocorrelation functions of quantities such as particle position and velocity of 
confined particles decay exponentially 0, [2^ . Additionally, the quantizations of classically chaotic billiards bear 
the signatures of classical chaos in, for example, the distributions of their energy levels, the "scarring" / "antiscarring" 
(increased/decreased wavefunction density that corresponds to classical unstable/stable periodic orbits, and the 
structure of nodal curves (at which the density vanishes) XL 2^ • 



Billiards with Mixed Dynamics 

Typical Hamiltonian systems exhibit mixed regular-chaotic dynamics (divided phase space). Examples of billiard 
geometries that yield such dynamics include non-concentric annuli (with circular boundaries) 29] and ovals [3|- For 
systems with mixed dynamics, however, it is difficult to exactly determine and describe the structure of the boundaries 
of stability islands. A container shaped like a mushroom provides an example of a billiard with divided phase space 
for which such precise mathematical analysis is feasible [J S • One reason it is natural to investigate mushroom 
billiards and their generalizations is that they can be designed in a precise manner so that they have an arbitrary 
number of integrable and chaotic components, each of which occupies the desired fraction of the phase space volume. 
Their study may thus lead to a better understanding of Hamiltonian systems with divided phase space. 

A circular mushroom billiard consists of a semicircular cap with a stem of some shape attached to the cap's 
base. Circular mushroom billiards provide a continuous transition between integrable (semi) circular billiards and 
chaotic (semi) stadium billiards. In circular mushrooms, trajectories that remain in the cap are integrable, whereas 
those that enter the stem are chaotic (except for a set of measure zero). The precise characterization of all such 
trajectories is completely understood, and one can change the dimensions of the mushroom to controllably alter the 
relative volume fractions in phase space of initial conditions leading to integrable and chaotic trajectories. Another 
interesting property of circular mushroom billiards, which was demonstrated recently, is that they exhibit "stickiness" 
in their chaotic trajectories (characterized by long tails in recurrence-time statistics) even though they do not possess 
hierarchies of K AM islands [3 • 

In this section, we investigate (axially symmetric) elliptical mushroom billiards, which generalize the elliptical 
stadium billiard 0? ^ind consist of a semiellipse with a stem extending from the center of its major (horizontal) 
axis. First, we show with an example plot (Fig.^ how to visualize the dynamics in both configuration space and phase 
space. Utilizing mushrooms with several stem geometries (rectangular, triangular, and trapezoidal), we subsequently 
examine symmetric and asymmetric periodic orbits and illustrate the birth and death of KAM islands that occur as 
the stem height is increased. As an example, we derive analytically the stem height at which the period-2 orbit along 
the mushroom's axis of symmetry becomes unstable. 

To study the dynamics of elliptical mushroom billiards, we conducted extensive numerical experiments using software 
that we wrote and have made publicly available [2^. To illustrate subsequent plots, we depict in Fig. ^ the billiard 
system in configuration and phase space for the case of a rectangular stem. We then show the dynamics of two of this 
billiard's trajectories in Fig. [2l The phase space of this billiard is divided, as it has both regular and chaotic regions. 
In mushroom billiards with rectangular stems, we observe stable periodic orbits of arbitrary length. 

We also study elliptical mushrooms with triangular (see Fig. (HI) and trapezoidal stems (no figures shown). When 
outside the "extended stem" of the mushroom, trajectories behave in the same manner regardless of the shape of the 



FIG. 1: (Color online) Configuration space (a) [in continuous time] and phase space (b) [in discrete time] plots for an elliptical 
mushroom billiard with a rectangular stem. The semimajor and semiminor axes of the mushroom's elliptical cap have respective 
lengths of 2 and 1. The stem has height 3 and width 1. In panel (b), the horizontal axis indicates the location of a boundary 
collision, as measured by the arclength in the clockwise direction from the right-most point of the table. The vertical axis 
indicates the incident angle of the particle's collision with the boundary. To illustrate how to read subsequent figures in this 
paper, we have plotted a trajectory (in black) that enters the mushroom's stem and another one (in red/gray) that does not. 
We label the regions in (b) to help identify their corresponding boundary arcs in (a). (Regions matched to different arcs are 
separated in (b) by vertical lines, which correspond to singular points of the billiard boundary.) For example, collisions against 
the elliptical arc of the mushroom's cap occur in region 6 (the right-most region) of phase space. 
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(a) « (b) 

FIG. 2: Periodic orbits of the elliptical mushroom from Fig.Q (a) Trajectory near a stable, symmetric periodic orbit of period 
34 (medium gray), (b) Trajectory near a stable, asymmetric periodic orbit of period 9 (light gray). 

stem. (The extended stem of a mushroom billiard consists of the stem itself plus the portion of the hat inside the 
largest geometrically similar semiellipse concentric to the hat that touches an upper corner of the stem.) Differences 
in dynamics for mushrooms with different stems manifest only for trajectories that enter the extended stem (and 
hence the stem per se). We remark that mushroom billiards with trapezoidal stems possess nontrivial periodic orbits 
that never leave the stem. In contrast, there are no such periodic orbits for mushroom billiards with triangular stems 
and only trivial orbits of this kind for mushrooms with rectangular stems. 

Symmetric mushroom billiards with rectangular, triangular, and trapezoidal stems all have stable symmetric and 
asymmetric periodic orbits that traverse both the hat and the stem. For triangular and trapezoidal stems, asymmetric 
periodic orbits that enter the stem are one of two types: (1) the trajectory collides against the stem near one of its 
bottom corners; or (2) the trajectory collides against a line segment that is perpendicular to the current link of the 
trajectory. (In the latter case, the particle then turns around and begins to trace its path backwards.) Rectangular 
stems, however, only have orbits of type 1. These orbits return to the cap from the stem at the same angle (with a 
vertical reflection) at which they left, with resulting boundary collisions in the cap as if the trajectory had bounced 
once off a linear segment perpendicular to the trajectory rather than bouncing several times in the stem (see, for 
example. Fig. Eb). Such trajectories enter the stem, bounce around against the stem's boundary, and then emerge 



FIG. 3: (Color online) Divided phase space and periodic orbits of an elliptical mushroom billiard with a stem shaped like an 
isosceles triangle. The semimajor and semiminor axes have lengths 2 and 1, respectively. The triangular stem has a height and 
width of 1. (a) Phase space showing multiple islands. The figure depicts 13 initial conditions, (b) Trajectory near a stable, 
symmetric periodic orbit of period 14 (shown in black), (c) Trajectory near a stable, asymmetric periodic orbit of period 10 
(shown in gray), (d) Phase space depicting the periodic orbits from (b) and (c). 

from it in the same fashion that they entered. By contrast, type 1 asymmetric periodic orbits in triangular stems do 
not have this property (see, for example, Fig.I^h). This difference in dynamics reveals the presence of periodic orbits 
in mushroom billiards with rectangular stems that are not present in ones with triangular stems. For symmetric 
periodic orbits, we observe a similar reflection property for all three types of stems, except that there is an additional 
horizontal reflection in the trajectory's return from the stem to the cap. Such trajectories (see Figs.|2ji andlHt)) result 
in periodic orbits that are present in all three families of mushroom billiards. No matter what they do inside the stem, 
such orbits leave the stem with the same angle (ignoring sign) with respect to the mushroom's axis of symmetry that 
they had when they entered. 

Numerous bifurcations occur with changes in stem height, as KAM islands rapidly proliferate and change con- 
figurations. For instance, using triangular stems of width 1 and heights ranging from 0.1 to 2.0 in increments of 
0.1, we observed a large variation in the number of central islands appearing vertically down the left side of phase 
space (corresponding to collisions against the billiard stem). For example, there are 6 such islands in Fig. [^h- In 
some configurations, the islands correspond to two different trajectories: an asymmetric trajectory and its horizontal 
reflection. In other configurations, one symmetric trajectory hits all the islands. We observed similar phenomena for 
trapezoidal stems with lower width 1.5, upper width 0.5, and heights ranging from 0.1 to 2.0 (in increments of 0.1). 

The bifurcation from stability to instability for the periodic vertical trajectory in symmetric elliptical mushrooms 
(for stems with a neutral arc at the bottom) occurs when the distance between the bottom of the stem and top of the 
cap is twice the radius of curvature of the cap. (More generally, this result applies to period-2 orbits which collide 
alternatively with a focusing elliptical arc and a straight (neutral) arc, which also occurs, for example, in the so-called 
"flattened" elliptical semistadium billiards with two foci H^.) The eigenvalues of the stability matrix (Jacobian) for 
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(a) ^ (b) 

FIG. 4: (Color online) Magnifications of phase space showing the disappearance of the central KAM island — which corresponds 
to the vertical period-2 trajectory passing through the axis of symmetry of the elliptical mushroom — as a result of varying the 
stem height. The example depicted has a semimajor axis of length 2, a semiminor axis of length 1, and a rectangular stem 
of width 1. The central island disappears at stem height 3 (when the trajectory becomes unstable). The figure panels show 
an enlargement of the portion of phase space corresponding to collisions with the elliptical arc. (a) Stem height 2.5. (b) Stem 
height 3. 

this period-2 orbit are 

Ai,2 = -i + |±^V-3-2z + z2, (1) 

which are complex conjugates when the stem height z < 3 and real (with one eigenvalue having magnitude greater 
than unity) when z > 3. This vertical trajectory thus becomes unstable at z = 3. This result is illustrated numerically 
in Fig. ^ which depicts an enlarged portion of phase space corresponding to collisions with the elliptical arc. The 
central island corresponding to the vertical trajectory shrinks as the stem height is increased until it finally disappears 
at height 3. 

Bunimovich proved that elliptical mushroom billiards with sufficiently long stems have one chaotic component 
and two integrable islands provided (1) the stem does not intersect the edges of the cap and (2) the stem does not 
contain the center of the cap's base Iq. (If the stem is not sufficiently long, other islands can also appear. One 
can thus observe such islands, for example, in elliptical stadium billiards jl^.) If condition (1) does not hold, the 
billiard contains one integrable island, formed by trajectories in the cap tangent to hyperbolae. If condition (2) does 
not hold, the elliptical mushroom contains only the island formed by trajectories in the cap tangent to ellipses. If 
both conditions are violated, then no integrable islands exist. In the limit of zero eccentricity, one obtains a circular 
mushroom billiard and the islands resulting from condition (2) disappear forever. Consequently, circular mushrooms 
contain no integrable islands if the stem intersects the edge of the cap. 

TWO-PARTICLE BILLIARDS 

We consider here the dynamics of two identical hard balls confined in containers whose one-particle billiard coun- 
terparts have integrable, chaotic, and mixed dynamics. In these systems, the confined particles collide elastically not 
only against the boundary but also against each other. We studied particles of various sizes to examine the "billiard" 
(low density) limit in which the particle radius vanishes. We find that although the two-particle dynamics are chaotic 
for all three situations, two-particle billiards with a chaotic or mixed billiard limit have different dynamics from those 
with an integrable billiard limit. 

The configuration space of a hard ball of radius r in a container Q C M'^ is equivalent to that of a point particle 
in a smaller container Qo C Q. The boundary dQo is formed by all the points q G int((5o) located a distance p from 
the boundary dQ. In the case of N noninteracting hard balls (that pass through each other and collide only against 
the boundary), the resulting configuration space is the direct product Qat := Qo x • • • x Qo times) and the phase 
space similarly consists of N copies of the one-particle phase space. 

For N interacting particles, it is difficult in general to explicitly describe the configuration space Qat, which can 
have a very complicated topology. For hard balls, the interaction potential is infinite if the particles collide and zero 
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FIG. 5: Distribution of incident angle for collisions between particles and boundary for two identical particles of varying radius 
confined in a circle of radius 1. The initial conditions for the three plots are shown in the last row of Tableland discussed 
in the text. We simulate this system through 200,000 total collisions of either particle with the boundary. (In this and all 
subsequent multiple-particle simulations, at most one data point is generated after each collision, as the incident angle remains 
the same if there are no particle-particle collisions.) (a) Particle radius 1/4. (b) Particle radius 1/10. (c) Particle radius 1/20. 
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TABLE I: Mean /i and standard deviation a of the incident angle for collisions between the particles and boundary for two 
particles of radius 1/4 confined in a circle of radius 1. The center of particle i has initial position (xi^yi) with initial velocity 
Vi at angle ai measured from the horizontal. (The choice of initial conditions is discussed in the text.) 

otherwise. To obtain the configuration space, one thus removes from Qn the N{N — l)/2 cylinders corresponding to 
pairwise interactions between particles 9]. (One obtains these cylinders by considering collisions between two fixed 
particles and allowing all others to move (without colliding) inside Q.) 

Two Particles in Geometries Corresponding to Integrable Billiards 

First, we consider two-particle dynamics in configurations which are integrable for point-particle billiards with a 
single particle. Using a confining circular boundary (of unit radius), we examined particles with radius 1/4, 1/10, and 
1/20 to study the dynamics as one approaches the underlying integrable billiard (i.e., as there are successively fewer 
collisions between the two particles). 

In Fig. 13 we show the distribution of incident angles for boundary collisions with two identical particles confined 
in a circle of radius 1. (The incident angle with the boundary is a natural observable for multiple-particle circular 
billiards, as it is a constant of motion for one-particle circular billiards. Because of the system's symmetry, the 
coordinate along the boundary is not essential.) We observe noticeable differences in the left tail of the distribution 
between the radius 1/4 particles and the two smaller particles, as tails in particle-particle versus particle-boundary 
collisions get heavier as the particle size becomes smaller. For the radius 1/4 particles, this distribution has a standard 
deviation of a ~ 0.3655 By contrast, it is about 0.4024 for radius 1/10 and about 0.4112 for radius 1/20. We also 
note that the collisions have a non-zero mean angle (about 0.5), as indicated in Tabled 

Although we obtain the same distribution of collision angles for all three particle sizes, we show in Fig. |H1 that one 
can separately track the two particles provided they are sufficiently large. As the particle radii are increased relative 
to the container radius, their motion becomes further constrained until they can always be distinguished (tracked) 
even in this ergodic setting. This is not true for general two-particle systems. In Fig. IHl we plot the incident angle 
of particle 1 versus the incident angle of particle 2 and color-code the figure based on their speed ratio Vi/v2. Bluer 
dots signify larger Vi/v2 > 1, greener dots signify smaller Vi/v2 < 1, and black dots signify Vi/v2 = 1. The colors 
are scaled from blue to green, so that brighter colors signify a larger speed disparity. For example, blue-black dots 
indicate that particle 1 is only slightly faster than particle 2. 

We choose the initial conditions for these simulations as follows: The center of particle 1 is selected uniformly at 



FIG. 6: Incident angle of particle 1 versus incident angle of particle 2 for data from Fig.|^ The colors signify the speed ratio of 
particle 1 to particle 2, with blue dots indicating vi > V2, green dots indicating vi < V2, and black dots indicating vi — V2- The 
colors are scaled so that collisions with speed ratios close to 1 are shown as nearly black and those with larger speed disparities 
are brighter, (a) Particle radius 1/4. (b) Particle radius 1/10. (c) Particle radius 1/20. As shown in panel (a), the particles 
can be tracked individually when they are sufficiently large. 

random from all points inside the circle that are not within 1/4 of the boundary of the circle. The center of the second 
particle, which must have distance at least 1/2 from the first initial coordinate and distance at least 1/4 from the 
circle's boundary, is chosen uniformly at random from its allowed points. We choose the angles uniformly at random 
between — tt and tt. We then select the first speed uniformly at random from [0, 2] and choose the second speed to 
ensure that the sum of the squares of the speeds is 4. 

We measure the locations of collisions between two particles in terms of their distance from the center of the table. 
The smaller the particle, the farther away from the origin one expects to find these collisions. For radius 1/4 particles, 
the distance from the origin peaks near 0.48 and drops off approximately linearly on both sides of the peak until it 
reaches near distances 0.28 and 0.71. For radius 1/10 particles, the collision location distribution peaks at roughly 
0.78 and drops off sharply to the right; it reaches near 0.35 and 0.9. For radius 1/20 particle, collisions peak near 
distance 0.88 from the center and their distribution again drops off sharply to the right, reaching near 0.38 and 0.94. 

We investigate the distributions of the number of total particle-boundary collisions between successive inter-particle 
collisions, the average intervening times between which are 1.8320, 4.8933, and 10.3158 for particles of radius 1/4, 
1/10, and 1/20, respectively. As expected, this distribution drops off more sharply for smaller particles. For example, 
for radius 1/4 particles, more than 90% of inter-particle collisions occur before there are 6 consecutive collisions of a 
particle with the boundary. For particles of radius 1/10 and 1/20, the number of consecutive collisions at this 90th 
percentile are 11 and 22, respectively. The modal number of boundary collisions between inter-particle collisions is 1 
in each case, with respective occurrence probabilities of 0.3361, 0.2495, and 0.1586 for particles of radius 1/4, 1/10, 
and 1/20. 

Now consider a pair of particles, again confined in circular containers, with nearly the same and nearly opposite 
initial velocities, so that we are perturbing from periodic orbits in which the particles collide against each other along 
the same sequence of chords of the circle. We show plots of the incident collision angle of one particle versus the 
other for the former case in Fig.E^-d, with slightly different initial ^/-coordinate values in each plot, and for the latter 
case in Fig. Et. (Given the same initial angles, one obtains essentially the same plots if the other initial conditions 
are changed.) Observe the change in dynamics from panel (a) to panel (d). We show only one panel for the case of 
nearly opposite velocities because the dynamics do not change with different ^(0). The mean incident angle against 
the boundary (averaging over both particles) is nonzero in (a)-(d) but zero in (e). Observe that the region in (d) is 
convex and does not develop any "wings" like the ones in panels (a)-(c), so that the ranges of incidence angles of the 
particles are very small for sufficiently large offset ^(0). Using the same color-coding as in Fig. El we see in the case of 
particles with nearly opposite initial velocities that if the incident angle of particle 1 is greater than that of particle 
2, then particle 1 has a lower speed than particle 2, as the line indicating equal incident angles sharply divides the 
figure into four monochromatic quadrants (observe the symmetry between the two particles). This is demonstrated 
in panel (f), in which all collisions with vi/v2 > 1 are shown in blue (dark gray) and all collisions with V2/V1 < 1 
are shown in green (light gray). Panel (g) shows that in this case, the magnitude of the angle ratio between the two 
particles is approximately equal to that of the velocity ratio, so that one needs less information to effectively describe 
the state of the system than would otherwise be necessary. 
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FIG. 7: Two particles of radius 1/10 in a circular billiard of unit radius starting in nearly the same or nearly opposite directions. 
For each initial condition except that in panel (d), we numerically simulate the collisions between the two particles for 50000 
total collisions. (In all of our simulations, we use the terminology 'total collisions' to indicate all the collisions between two 
particles and between either particle and the boundary.) We run our simulation through 500000 such collisions in (d) to confirm 
the accuracy of the result. Each point represents one collision between the two particles. The horizontal axes depict the incident 
angle of particle 1, and the vertical axes depict the incident angle of particle 2. Both particles have an initial speed of 1. The 
first particle in all panels has initial horizontal coordinate xi(0) = 0.3, and the second particle always has initial horizontal 
coordinate X2(0) = —0.3. In all cases, both particles have the same initial vertical coordinate y{0). The colors are defined as 
in Fig. IHl Panels (a)-(d) depict the results for particles with nearly the same initial angles [^i(O) = 0, ^2(0) = 0.1], and panel 
(e) depicts the result for particles with nearly opposite initial angles [^i(O) = tt, ^2(0) = 0.1]. (We obtain the same results for 
smaller initial angle differences.) For the case of nearly opposite angles, the dynamics are the same for different y{0), so we 
only show one plot. Panels (a) and (e) have ^(0) = 0.2, (b) has y{0) = 0.4, (c) has y{0) = 0.6, and (d) has y{0) = 0.8. (Color 
online) Panel (f) shows the same simulation as panel (e), except all collisions with vi > V2 are shown in blue (dark gray) and 
all collisions with vi < V2 are shown in green (light gray). Panel (g), also for this simulation, shows a plot of the magnitude 
of the angle ratio versus the magnitude of the velocity ratio, indicating that they are approximately equal for particles with 
almost opposite initial angles. 



Two Particles in Geometries Corresponding to Chaotic Billiards 



Diamond billiards are dispersing with boundaries consisting of arcs of four circles which cross at the vertices of a 
square. We studied two-particle diamond billiards both with and without tangencies between the circular arcs. In 
Fig.lHtb, we show the distributions of their incident angles for particle-boundary collisions. They differ markedly from 
those observed for two hard balls in circular containers (see Fig. [5|). In particular, the mean incident collision angle 
is for two-particle diamond billiards but nonzero for two-particle circular billiards. In the diamond containers, the 
distributions fall off from the central peak to at about ±7r/2. The distributions for two-particle diamond billiards are 
also wider than they are for circular ones. Also unlike the case of circles, we observe here a higher density of collisions 
near the corners, so that the Poincare map includes a clustering of points near the vertical lines corresponding to 
particle collisions with these corners. 



Two Particles in Geometries Corresponding to Billiards with Mixed Dynamics 



To examine two-particle billiards whose one-particle counterparts have mixed dynamics, we again consider containers 
shaped like mushrooms. In Fig. [St, we show the distribution of incident angles for particle-boundary collisions in the 
case of a circular mushroom. This distribution is reminiscent of that for the diamond geometries in Fig.jStb. In these 
two-particle simulations, we examined initial conditions with both particles in integrable regions of the corresponding 
one-particle billiard, both particles in chaotic regions, and one particle in each of these types of regions. In all cases, 
we obtain the same results as with the diamond geometry (which, again, has a fully chaotic billiard limit). We also 
observe these same dynamics for elliptical mushrooms with very eccentric caps. (We considered examples with the 
stem attached to the foci, so that only one family of caustics is present in the one-particle billiard.) If islands are 




10 



[a) [D) [c) 

FIG. 8: Distribution of the incident angle of (several thousand) particle-boundary collisions for two hard balls of radius 1/4 
confined in (a,b) diamond containers composed of arcs of radius- 1 circles and (c) a circular mushroom container. The diamond 
in (a) has tangencies, as it consists of four quarter- circle arcs. The diamond in (b) does not have tangencies, as it consists of 
four arcs that each constitute l/8th of a circle. 







(a) (b) (c) ~ (d) 



FIG. 9: Distribution of incident angle for particle-boundary collisions for three identical particles of radius .0816 confined in a 
circle of radius 1. In each case (covering 100000 total collisions), we initialize all three particles with unit speed and the same 
y coordinate (which is varied in the simulations). For each simulation, the initial horizontal coordinates of the three particles 
are xi(0) = 0.3, X2(0) = 0, and X3(0) = -0.3. The initial angles are 6*1(0) = 6*3(0) = and 6*2(0) = 0.1. We list the observed 
means and standard deviations in Table UTI The initial vertical coordinates are (a) ^(0) = 0, (b) y{0) = 0.2, (c) y{0) = 0.4, and 
(d) y(0) = 0.6. 

present, they are probably very small in size. 

Three- and Four-Particle Billiards 

Because KAM tori divide phase space in two dimensions but not in higher dimensions, it is crucial to consider 
containers that confine three or more particles. One may thus expect some differences in the dynamics when one 
considers more than two particles. 

To see that signatures of integr ability remain even with additional particles, we considered three hard disks of 
radius 0.0816 in a circle of unit radius. (The area occupied by the three particles in these simulations is almost the 
same as that occupied by the radius 1/10 particles in the corresponding two-particle simulations.) In particular, we 
investigated initial conditions in which all three particles are collinear and initially traverse almost the same chord 
of the circle (i.e., their initial angles are almost the same). As in the two-particle case, one of the balls has a small 
angular displacement (of 0.1) with respect to the others and the common ^/-coordinate is varied exactly as for the 
two-particle simulations. As shown in Fig. El we again observe a non-zero mean in the distribution of incident angles 
for particle-boundary collisions. (For a given simulation, the distribution for each of the individual particles is roughly 
the same, so we show the distributions for all three particles collectively.) 

We also observe a non-zero mean in the distribution of incident angles in a similar numerical experiment with four 
hard balls in a circle. Our observations are summarized in Table ITTI 



CONCLUSIONS 



In this paper, we examined the dynamics of one-particle and few-particle billiards. These situations can be compared 
in the context of noninteracting versus interacting hard balls in a given container. (The "billiard limit" discussed 
in the paper in which the particle size becomes smaller also provides a means for comparison.) For example, two 
noninteracting hard balls in a circular container are integr able, as this system just consists of two integr able billiards. 
We showed using numerical simulations that two interacting hard balls in a circular container bear signatures of 
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y(o) 


^l2 


M3 


/i4 


cr2 


crs 







0.0083 


0.0088 


-0.0085 


0.4256 


0.4863 


0.5195 


0.2 


0.2143 


0.1893 


0.1894 


0.4788 


0.4899 


0.5189 


0.4 


0.4829 


0.4529 


0.4443 


0.4283 


0.4627 


0.4793 


0.6 


0.8241 


0.7957 


0.7852 


0.2485 


0.3132 


0.3490 



TABLE II: Mean jij and standard deviation Gj for the incident angle of particle-boundary collisions for j — 2,3,4 hard balls 
in a circular container. In all simulations (of 100000 total collisions), each particle has initial speed 1 and initial y-coordinate 
^(0). The other initial conditions for j — 2 [j — 3) are shown in Fig. [7| (Fig. 1^. For j — 4, we considered circular particles 
of radius 0.0707. For each simulation, the initial x-coordinates are xi(0) = 0.3, X2(0) = 0.1, X3(0) = —0.1, and xa{0) — —0.3. 
The initial angles are ^i(O) = ^3(0) = ^4(0) = and ^2(0) = 0.1. Observe that the mean decreases and the standard deviation 
increases with the additional particles. 



integrability even though the system is chaotic and that its dynamics differs from that of two particles in a diamond 
container (which is chaotic even in the one-particle case). Also, two noninteracting particles in a mushroom show 
mixed dynamics, whereas two interacting particles in a mushroom container seem to show the same dynamics as two 
interacting particles in a diamond. 

Consequently, we have shown that the shape of a container confining a system of interacting particles does matter, 
despite the fact that this aspect of the dynamics in the behavior of a gas of particles seems to be essentially neglected 
in investigations of these systems. In fact, the container shape may be of comparable importance to that of the 
interaction potential. Although the importance of the shape of the container (the so-called "billiard table") for 
systems of noninteracting particles is universally acknowledged, a "silent consensus" still exists that if the interaction 
potential is nontrivial, then a system of interacting particles must be ergodic in the thermodynamic limit. We have 
shown in this paper, however, that this is not true for systems of a few interacting particles, which have thus far 
received little attention. With the emergence of nanosystems and new experimental techniques, such systems are now 
among the most important ones for various applications 0, . 

Exactly solvable billiard systems demonstrating all three possible types of behavior (integrability, chaos, and mixed 
dynamics) motivate natural container shapes to study the dynamics of a few nontrivially interacting particles. While 
the first small step has been made in this paper, most of the problems in this area remain completely open; we mention 
a few of them in passing. First, it is well-known in statistical physics that the high-density limit is more regular than 
the low-density limit (where, for example, logarithmic terms appear in the asymptotics of transport coefficients |l^). 
We have observed this in our numerical studies as well. In fact, it seems to be a general phenomenon for which a 
theory should be developed without the ubiquitous assumption that the number of particles is a large parameter. 
As usual, the temporal asymptotics of both correlation functions and Poincare recurrences should be investigated in 
various situations, including all three types of behavior in billiards (noninteracting particles) and for broad classes 
of potentials. One of the key questions is what happens to the relative volume of KAM islands when the number of 
particles tends to infinity [sf. Again the "silent consensus" (i.e., a belief without proofs or convincing demonstrations) 
is that this relative volume tends to zero. However, it has been demonstrated that this is at least not always true 
Hence, it is not at all unlikely that nonlinear dynamics will reveal some new surprises in further investigations of this 
phenomenon. 



Acknowledgements 

SL and MAP acknowledge support provided by an NSF VIGRE grant awarded to the School of Mathematics at 
Georgia Tech. The work of LB and SL was partially supported by NSF grant DMS 0140165. LB also acknowledges 
the support of the Humboldt Foundation, SL acknowledges funding from a Georgia Tech PURA award, and MAP 
acknowledges support from the Gordon and Betty Moore Foundation through Caltech's Center for the Physics of 
Information. We also thank Nir Davidson for useful conversations and the anonymous referee for suggestions that led 
to significant improvements in the paper. 



[1] B. J. Alder and T. E. Wainwright. Decay of velocity autocorrelation function. Physical Review A, 1:18-21, January 1970. 
[2] E. G. Altmann, A. E. Motter, and H. Kantz. Stickiness in mushroom billiards. Chaos, 15(3):033105, September 2005. 



12 



[3] M. F. Anderson, A. Kaplan, T. Griinzweig, and N. Davidson. Decay of quantum correlations in atom optics billiards with 

chaotic and mixed dynamics. arXiv: quant-ph/0404118 October 2004. 
[4] D. V. Anosov and Ya. G. Sinai. Some smooth ergodic systems. Russian Mathematical Surveys, 22:103-167, 1967. 
[5] V. I. Arnold. Proof of A. N. Kolmogorov's theorem on the preservation of quasiperiodic motions under small perturbations 

of the Hamiltonian. Russian Mathematical Surveys, 18(5):9-36, 1963. 
[6] V. I. Arnold. Small divisor problems in classical and celestial mechanics. Russian Mathematical Surveys, 18(6):85-192, 

1963. 

[7] M. V. Berry. Regularity and chaos in classical mechanics, illustrated by three deformations of a circular billiard. European 

Journal of Physics, 2:91-102, 1981. 
[8] L. A. Bunimovich. Mushrooms and other billiards with divided phase space. Chaos, ll(4):802-808, December 2001. 
[9] L. A. Bunimovich. Kinematics, equilibrium, and shape in Hamiltonian systems: The 'LAB' effect. Chaos, 13(3):903-912, 

2003. 

[10] P. Cvitanovic, R. Artuso, R. Mainieri, G. Tanner, and G. Vattay. Chaos: Classical and Quantum. Niels Bohr Institute, 

Copenhagen, 11 edition, February 2005. ChaosBook.org. 
[11] S. De Bievre, P. E. Parris, and A. Silvius. Chaotic dynamics of a free particle interacting linearly with a harmonic oscillator. 

Physica D, 208(l-2):96-114, August 2005. 
[12] M. De Leener and M. Resibois. Classical Kinetic Theory of Fluids. John Wiley and Sons, New York, NY, 1977. 
[13] G. Del Magno and R. Markarian. Bernoulli elliptical stadia. Communications in Mathematical Physics, 233(2) :21 1-230, 

February 2003. 

[14] Ch. Dellago, H. A. Posch, and W. G. Hoover. Lyapunov instability in a system of hard disks in equilibrium and nonequi- 

librium steady states. Physical Review E, 53(2): 1485-1501, February 1996. 
[15] N. Friedman, A. Kaplan, D. Carasso, and N. Davidson. Observation of chaotic and regular dynamics in atom-optics 

billiards. Physical Review Letters, 86(8):1518-1521, February 2001. 
[16] P. Garrido and G. Gallavotti. Billiards correlation functions. Journal of Statistical Physics, 76:549-586, 1994. 
[17] M. C. Gutzwiller. Chaos in Classical and Quantum Mechanics. Number 1 in Interdisciplinary Applied Mathematics. 

Springer- Verlag, New York, NY, 1990. 
[18] D. M. Harber, H. J. Lewandowski, J. M. McGuirk, and E. A. Cornell. Effect of cold collisions on spin coherence and 

resonance shifts in a magnetically trapped ultracold gas. Physical Review A, 66:053616, 2002. 
[19] A. Katok and B. Hasselblatt. Introduction to the Modern Theory of Dynamical Systems. Cambridge University Press, New 

York, NY, 1995. 

[20] A. N. Kolmogorov. On conservation of conditionally periodic motions under small perturbations of the Hamiltonian. Dokl. 

Akad. Nauk. SSSR, 98:527-530, 1954. 
[21] A. Kudrolli, V. Kidambi, and S. Sridhar. Experimantal studies of chaos and localization in quantum wave functions. 

Physical Review Letters, 75(5):822-825, July 1995. 
[22] S. Lansel and M. A. Porter. A GUI billiard simulator for Matlab. nlin.CD/0405033, May 2004. 

[23] V. Lopac, I. Mrkonjic, N. Pavin, and D. Radic. Chaotic dynamics of the elliptical stadium billiard in the full parameter 

space. nlin.CD/0507014, July 2005. 
[24] S. W. MacDonald and A. N. Kaufman. Spectrum and eigenfunctions for a Hamiltonian with stochastic trajectories. 

Physical Review Letters, 42(18):1189-1191, April 1979. 
[25] S. W. MacDonald and A. N. Kaufman. Wave chaos in the stadium: Statistical properties of short-wave solutions of the 

Helmholtz equation. Physical Review A, 37:3067, 1988. 
[26] C. M. Marcus, A. J. Rimberg, R. M. Westervelt, P. F. Hopkins, and A. C. Gossard. Conductance fluctuations and chaotic 

scattering in ballistic microstructures. Physical Review Letters, 69(3):506-509, July 1992. 
[27] J. Moser. On invariant curves of area-preserving mappings of an annulus. Nachr. Akad. Wiss. Cottingen Math. Phys. Kl., 

2:1-20, 1962. 

[28] H. A. Posch and R. Hirschl. Simulations of billiards and of hard body fluids. In D. Szasz, editor. Hard Ball Systems and 
the Lorentz Cas, volume 101 of Encyclopaedia of Mathematical Sciences, pages 279-314. Springer- Verlag, Berlin, Germany, 
2000. 

[29] N. Saito, H. Hirooka, J. Ford, F. Vivaldi, and G. H. Walker. Numerical study of billiard motion in an annulus bounded 

by non- concentric circles. Physica D, 5:273-286, 1982. 
[30] Ya. G. Sinai. On the foundation of the ergodic hypothesis for a dynamical system of statistical mechanics. Dokl. Akad. 

Nauk. SSSR, 153:1261-1264, 1963. 
[31] Ya. G. Sinai. WHAT IS. ..a billiard. Notices of the American Mathematical Society, 51(4):412-413, April 2004. 
[32] Ya. G. Sinai. Dynamical systems with elastic reflections. Russian Mathematical Surveys, 25(N2): 137-188, 1970. 
[33] S. Smale. Differentiable dynamical systems. Bulletin of the American Mathematical Society, 73:747-817, 1967. 
[34] J. Stein and H. J. Stockmann. Experimental-determination of billiard wave- functions. Physical Review Letters, 68(19):2867- 

2870, 1992. 

[35] H. J. Stockmann and J. Stein. Quantum chaos in billiards studied by microwave-absorption. Physical Review Letters, 
64(19):2215-2218, 1990. 

[36] C. Zhang, J. Liu, M. G. Raizen, and Q. Niu. Quantum chaos of Bogoliubov waves for a Bose-Einstein condensate in 
stadium billiards. Physical Review Letters, 93:074101, 2004. 



